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Abstract 

An asymptotic solution of the Boltzmann equation is developed for 
ICR absorption, without restrictions on the ion-neutral collision fre- 
quency or mass ratio. Velocity dependence of the collision frequency 
causes deviations from Lorentzian line shape. 
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1. Introduction 
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Most of the work on the theory of ion cyclotron resonance (ICR) 
absorption spectra has been restricted to the case of velocity-inde- 
pendent ion-neutral collision frequencies (Maxwell model), or has in- 
volved the essentially equivalent assumption that the average of a 
product is equal to"t:he product of the averages [1-6], These theories 
all yield Lorentzian line shapes for the steady-state absorption in a 
weakly ionized gas , provided that the effects of walls and chemical 
reactions can be ignored, and that the collision frequency is much 
less than the cyclotron frequency. The Maxwell -model restriction has 
been removed for electron absorption, but only at the cost of the new 
restriction of ion mass much less than neutral mass (Lorentz model) 
[7,8]. For this model it was found that the line shape depended 
strongly on the velocity dependence of the collision frequency. The 
problem could also be solved for the special case of ion mass much 
greater than neutral mass (Rayleigh model) [9], but, as shown below, 
the line shape will again be Lorentzian. 

The purpose of this communication is to lift the restrictions on 
ion-neutral collision frequency and mass ratio, so that deviations from 
Lorentzian line shape can be interpreted in terms of the nature of the 
ion-neutral force law. 

2, Theory and Results 

The fundamental assumption is that the ion density , n , is much 
less than the density of neutral gas molecules, N , so that ion-ion 
collisions can be neglected compared to ion-neutral collisions. The 
collision term of the Boltzmann equation then becomes linear, and can 
be expanded in a series and the resulting equtions solved by a moment 
method [10,11], without restrictions on collision frequency or mass 
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ratio. In order illustrate the dependence of the ICR line shape on 
the ion-neutral force law, we consider a very simple situation. Ion- 
neutral collisions are assumed elastic. A uniform static magnetic field 
B points in the y-direction, and a uniform alternating electric field 

■W\ 

E = E sin wt points in the z-direction. The ion density is assumed 
-v* /**>o 

uniform so that boundaries are neglected, and the ion-neutral colli- 
sion frequency is assumed small compared to the cyclotron frequency 
to^ = qB/rn , where q and m are the charge and mass of the ion. 
Finally, we seek only a steady-state solution, so that we can neglect 
the initial conditions of the ions, which decay exponentially in time. 
The theory is not limited to this specialized situation, and the fore- 
going restrictions can be relaxed (except for n << N and elastic 
collisions), although at the price of some additional mathematical 
complication. 

Solutions of the linearized Boltzmann equation can be obtained 
by expanding the collision operator in terms of a complete set of 
orthogonal functions and forming moments [10]. We use the Burnett 
functions, the eigenfunctions of the collision operator for the Maxwell 
model, 

■ •£’ ■ = £ O" 2 ' •> ' 1 "‘ • <“ 


c = v (m/2kT)^^^ , cos 0 = c /c , 


( 2 ) 


(r) 


where a Son ^ ne polynomial, Pj> m ’ an associated Legendre 

polynomial, v^ the ion velocity, and T the neutral gas temperature. 
The infinite set of moment equations must be solved by some truncation 
and iteration procedure; using the scheme in [11,12] we obtain for the 
p-th approximation the expression. 
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where the pointed brackets signify averages over the ion distribution 

function. The coefficients b (£m) are matrix elements of the lin- 

rs 

earized Boltzmann collision operator for multicomponent gas mixtures with 

(r ) 

respect to the , and are defined and tabulated in Ref. [12-14], 

Dependence on composition of the mixture is entirely in the b (£m) , 
which are linear in the mole fractions. The summation in the last term 


begins at s^ = r-p+1 or 0 , whichever is larger, and includes only 
those matrix elements less than p units off the diagonal. 


The iteration begins at p=l 


mat ion and solving for 



by neglecting the terms in the sum- 
, and then proceedings to p=2,3,etc. 


The mean power absorption A' from the electric field is the time aver- 


age of 


A(t) = nq^v(t) y » £(t) = nqE^ (v^Ct) ) sin art . (4) 

Thus we must solve eq. (3) for\i^^(t)/ = /c (t)\ . For p=l 

'10 / p \ z / p 

we obtain two coupled differential equations from eq. (3), and for 
p=2 we obtain eight partially coupled differential equations. Solu- 
tion of these increasingly large sets of differential equations leads 

— 2 

to a series for A in powers of (E/N) . The first iteration neglects 

all off-diagonal matrix elements and yields the usual Lorentzian line 
shape for A [1,2], This result will be exact for very weak electric 
fields, or for all electric fields if the off-diagonal matrix elements 
are zero. This is the case for the Maxwell model and the Rayleigh model, 



which thus always give Lorentzian line shapes. 

Deviations from Lorentzian line shape first appear in the second 
iteration, which yields 


A = 


2 r 2 

f I 

1 


nq E 
1 o 




2 


4mv 

U + n J 



1 + 


C 1 C 2 & 


(1 + ft 2 ) (c 2 + n 2 ) 


+ 0 (£ 7 ) + ... 


(5) 


where v = Nb^dO) is the zero-field mean collision frequency of the 
ions for momentum transfer to the neutrals in the mixture, and 


fl = (w - w )/v , 

(6) 

f = (qE o /mv) (m/2kT) 1/2 , 

(7) 

. 5b oi (1 °) . b 01 (10) 

C 1 " 3b u (00) + 6b Q0 (20) ’ 

(8) 

C 2 = b 11 ^ 1 °)/ b oo (l0) * 

(9) 


A physically interesting feature of the above formula is the occurrence 
of b^ClO) , which depends strongly on the nature of the ion-neutral 
force law. For ions in a single gas of mass M the expressions for 
c and c^ become 

1 A A J* n 

c 1 = - (6C“ - 5) (M/m) (11m + 6MA~) (5m + 3MA") , (10) 

c_ ^ (3m 2 + M 2 + -f- mMA ) (m + M) 2 , (11) 

/ o 

A A 

where C and A are conventionally defined [15] ratios of colli- 
sion integrals. The crucial part coming from b (10) is the factor 

j’{ # 

6C - 5 , which also plays a crucial role m the theory of thermal dif- 
fusion in gases. Discussions of the relation of thermal diffusion to 
molecular collisions and force laws can therefore be taken over directly 
and applied to the relation between ion-neutral interactions and devia- 
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tions from Lorentzian line shape [16,17], 

As an illustration of the effects on line shape of a repulsive ion- 
neutral force law, we consider the case of rigid spheres and equal 
masses (m=M) , for which c^ = - 17/96 and c^ = 7/5 . The results 
are shown in fig. 1 for several values of , The line is Lorentzian 
for £ = 0 , but the absorption is less for £ > 0 because the effec- 
tive collision frequency increases with increasing ion energy for repul- 

-4 

sive potentials "harder" than r . The effect is greatest at the 
line center because the ion energy is greatest at resonance. These re- 
sults are qualitatively the same for any gas temperature. 

For a more realistic ion-neutral interaction, such as a long-range 
attraction plus a short-range repulsion, the effective collision fre- 
quency has a more complicated dependence on field strength and gas tem- 
perature. For low gas temperatures the collision frequency may first 
decrease with increasing field strength, pass through a minimum, and 
then increase. For high gas temperatures it may increase from the be- 
ginning. The effect on the line shape is qualitatively illustrated in 
fig. 2. The qualitative effect can be inferred from the above consi- 
derations without explicit computations, just as the static field de- 
pendence of the ion mobility at a given gas temperature can be inferred 
from the temperature dependence of the zero- field mobility, as first dis- 
cussed by Kihara [10,11,14]. 

The foregoing remarks can be made mathematically explicit by ex- 
pressing the absorption line shape in terms of the ion mobility K in 
a static electric field. Eliminating £, between eq. (5) and the cor- 
responding expression for K to the same order of approximation [11,12,14], 
we obtain 

nqE 2 K(0) 

A = ^ i -f 

4(1 + ft ) 


2 

C 2 °3 

24(1 + ft 2 ) (c 2 + ft 2 ) 


» 


( 12 ) 
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where K(0) is the zero-field limit of the mobility, so that 
fl = (w - o> c ) q/mK(0) , and 

c 3 = 3 |t> 11 (°0) + ^ 0b oo (20) ]/[ 4b ll (00) + 5b 00 (2 °^j » (13) 

s V 

which for single gases reduces to = (11m + 6MA )/(3m + 4MA ) . 

Both c 2 and c 3 are virtually independent of the ion-neutral force 
law, which enters only weakly through A . As written, eqs. (12) and 
(13) are valid for multicomponent mixtures of neutral gases. 

3. Discussion 

Since ICR lines are often used as a probe for chemical analysis 

[18], it is important to note that A remains proportional to n 

even when the line shape is non-Lorentzian. 

— 2 

Since the present results give A as a function of ft , the 
line shape is predicted to be symmetric about 8=0. Skewness of 
line shape [19] therefore cannot be attributed to the nature of the 
ion-neutral force law. A similar remark applies to the central dips 
sometimes seen in ICR lines. 

The convergence of the expansion in may be slow, and eq. (12) 

is likely to have a much larger range of usefulness than eq. (5), if 
one judges by the analogous expressions for static-field ion mobility 
and diffusion [12]. 

The fundamental limitations of the present theory are n << N 
and elastic collisions. The other restrictions are readily removable 
in principle, the cost being only mathematical complexity. 
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Figure Captions 

Fig. 1. Reduced ICR absorption A(ft,£) as a function of ft = (to - w c )/ v 
for several values of reduced electric field strength £✓ , for the model 
of rigid spheres of equal mass (m = M) . The curve is Lorentzian for 

£= o . 

Fig. 2. Qualitative behavior of A (ft,£) for a potential consisting of 
a rigid sphere with long-range attraction. The plots are made for 
£ = 1 at several gas temperatures. The curve for T = 00 is the same 
as the curve for £ = 1 in fig. 1. 
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